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Abstract 

Maslov indices in periodic-orbit theory are investigated using phase 
space path integral. Based on the observation that the Maslov index is 
the multi- valued function of the monodromy matrix, we introduce a gener- 
alized monodromy matrix in the universal covering space of the symplectic 
group and show that this index is uniquely determined in this space. The 
stability of the orbit is shown to determine the parity of the index, and a 
formula for the index of the n-repetition of the orbit is derived. 



PACS numbers: 05.45.Mt, 02.40.Ma 

Over the past few decades a large number of studies have been made on 
periodic-orbit theory However, the Maslov index, which is an additional 
phase factor appearing in periodic orbit theory, doesn't seem to be thoroughly 
understood. For hyperbolic orbits, Robbins Q showed that these indices are 
the winding numbers which are defined by the invariant Lagrangian manifolds 
around these orbits. Moreover, it was conjectured that the same argument 
could be extended to elliptic orbits and more general orbits which have mixed 
stability. However, this is not the case since general orbits don't necessarily 
have such invariant manifolds around them. Furthermore, Brack et al. Q 
investigated periodic orbits in anisotropic harmonic oscillator (these orbits are 
elliptic), and showed the Maslov index of the n-repetition of the orbit /i„ is 
not equal to rt/ii. This result contradicts Robbins's conjecture, which leads to 

= "•Ml- 

In this letter, we propose a new approach to the problem which can be ap- 
plied to all periodic orbits, irrespective of the type of the stability. Our method 
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is based on phase space path integral of the partition function. The advantage 
of using phase space path integral is that it respects symmetries of the sys- 
tem: time translation symmetry and canonical invariance. In addition to global 
canonical transformations, our method allows us to use local canonical transfor- 
mations which act differently at each point of the orbit. They are regarded as 
gauge transformations with gauge group Sp(2rt, i?), and the topological nature 
of this gauge group plays essential role in our theory. Derivation of semiclassical 
trace formula using phase space path integral was given by Sakita et al.0, and 
similar analysis was made by Kuratsuji Q for coherent state path integral, but 
Maslov indices are neglected in these works. Maslov indices of semicalassical 
propagator (but not the periodic orbit) was investigated by Levit et al. ^, ^ 
using phase space path integral, and they showed that this index is given by 
the excess of the negative eigenvalues of second variation of action around the 
orbit. We treat Maslov indices of periodic orbits in the same way as theirs. 
There are two types of Maslov indices of periodic orbits: indices for the parti- 
tion function (time region), and that for the density of states (energy region). 
We calculate Maslov indices for the partition function first. Then it is easy to 
calculate Maslov indices for the density of states by taking into account the ex- 
tra phase contributions that arises in performing the Fourier-Laplace transform 
by the stationary phase approximation. 

The main idea of this paper is that the linearized symplectic flow around the 
orbit directly determines the Maslov index. Therefore we don't use Lagrangian 
manifolds in this paper. The linearized symplectic flow is represented by V{t) (a 
symplectic matrix with one parameter i), and V{t) is considered to be a curve 
on the group manifold of Sp(2n, R). We can see various properties of Maslov 
indices by continuous deformation of this curve. 

Our final results can be summarized in the following formula eq.(|l]) and 
eq.(H). The Maslov index of the orbit is 



where p is the numbers of elliptic pairs of eigenvalues, q is the number of inverse 
hyperbolic pairs of eigenvalues, and k is the winding number which will be 
defined later on. We can see that parity of the index is determined by p and q 
from this formula. The Maslov index of the n-repetition of this orbit is 



where Ui is the stability angle of the i-th elliptic eigenvalues. Note that /l(„ is 
not equal to n/zi if the orbit have elliptic eigenvalues. 

Let us start with the following identity for the density of states; 



^1 = p + q + 2k, 



(1) 





P{E) 



— Im g{E + ie), 



(3) 



TT 
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g(E) = Tr ^. 

E-H 

g{E) is obtained by Fourier-Laplace transform of the partition function: 

1 r°° 

g{E) = — dTe'^^/''Z{T). 

JO 

We use the phase space path integral of the partition function: 



Z(T) ^ Trexp 



-HT 

n 



VpDq exp 



- (p {pdq - Hdt) 



(4) 



(5) 



(6) 



Strictly speaking, we should write in a discrete form to obtain well-defined 
continuum limit. However, to make the discussion simple and show our essential 
idea, we use continuum notation throughout this paper. The detailed analysis 
of (H), including subtle problems concerning operator ordering and changing 
variables, will be reported in the forthcoming paper. [l3| . 

Let us evaluate by the stationary phase approximation. The stationary 
phase condition reads 



<j>{pdq — Hdt) 



0, 



(7) 



which leads to the Hamilton equations of motion. It follows from the periodic 
boundary condition of the path integral that the solutions of (|^) are periodic 
orbits. Thus the semi-classical approximation ZsdT) to Z{T) is given by a sum 
over the periodic orbits. 



Zsc{T)=Y,Kexp 



-R 



(8) 



Here B = §{pdq — Hdt) is the classical action of the periodic-orbit, and K is 
the contribution of the quadratic fluctuation around the orbit, 



A-^/l,xexp[.^>«|x(„l] 



x{t) = —^{Spi,Sqi,6p2, .. 
The explicit form of the second variation 6^R is 



(9) 
(10) 



dt x^(t) 

where J is a matrix of the form: 



6'R[x{t)] = - 



-J-~H {p{t),q{t))]x{t), (11) 



/ 



J 



(12) 



3 J 
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and j is a 2 X 2 matrix: 



J = 



-1 

1 



We rewrite (nTI) as 



5^R = / dt x^{t )Ip{t)x{t 



where t = t/T, and Ip is defined as 



^[t) = JD = j[--A{t) 



A{t)=TJH {p{t/T),q{t/T)). 



(13) 



(14) 



(15) 



(16) 

A is the connection for the displacement vector x, and D — d/dt — A is the 
covariant derivative. The parallel transport of x is defined by the equation 
Dx = (Figj^), which is the Hamiltonian equation of motion for x: 

jX^JH"{p{t),q{t))x. (17) 

Hereafter we slightly generalize the problem, and A is considered to be a gen- 
erator of symplectic group which satisfies A^ J + J A ~ 0. Then A is written 
as 

A{t) = JH{t), (18) 



where H is & symmetric matrix of even demension. (16) is the special case of 
(181). We deal with the quadratic path integral of the following form instead of 

(i 



Maslov index: 



K = I Vx cxp 



1 



— X 
2 



1 



, (19) 

^Det ilp 

Here Ip = J{d/dt - A), and A is defined in (|l|). The phase of ^ defines the 

(20) 



V^Det ilp ^/\DctIp\ 



M = 



(21) 



Here /^+(/i_) denote the number of the positive (negative) eigenvalues of Ip . 
In the continuum limit, both /i+ and become infinity, and the formula ( pl[ ) 
cannot directly define the Maslov index. However, we can evaluate the relative 
change of the index by deforming the operator Ip continuously and analyzing 
the points where the sign of the eigenvalue changes. (This is similar to the 
discussion given m |, I for the semi-classical propagator.) 

In the following, we calculate |Det Ip \ and ^ in (pO|). The procedure consists 
of three steps. 
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1. Diagonalization of D and calculation of |Dct Ip \ = |Dct D\. 

2. Diagonalization of Ip using the basis given in the first step. 

3. Calculation of /i by continuous deformations of Ip . 

The first step is essentially the same as the discussions given in [Q and ||] . The 
second and the third steps are first given in this paper. 

The first step is the calculation of |Det Ip\. For this purpose, we first diago- 
nalize the covariant derivative D: 

Dx = (^^^ - X = ex. (22) 

This equation is easy to solve because the time derivative has diagonal form. 
We introduce a symplectic matrix V(i) which stands for the flow around the 
orbit: 

DV{t) = 0. (23) 

If we fix V{Q) — 1, D and V has one to one correspondence, and V{1) (= M) is 
the monodromy matrix of the orbit. Let be the eigenvector of M, and ( be 
the stability angle of the orbit, 

= e-<|^. (24) 

Then the solutions of (^2|) are written as follows: 

a;^,„(t) = eM^Cnt)V{t)^^, (25) 

Dxi^.,i = e^,„a;^,„, (26) 

e^,„ = (C + 2mr)i. (27) 

If e~"' is the eigenvalue of M, e*^ also becomes the eigenvalue of M Q, and the 
absolute value of the functional determinant of Ip can be evaluated as a product 
of eigenvalues (27): 

IDet^l = iDetDl = ]J |(27rn)2 - Cf | , (28) 



n.l 




(29) 



Here we replace the divergent constant Jli^i ^Trn by 1. In discrete formalism, 
this term corresponds to 

1 -n- 7rn 

— IT 2 sin — = 1. (30) 

n=l 
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For n ^ N, J^^j 2 sin 7rn/7V ~ J|^27rn/iV, therefore (|3^) is considered to be a 
normalized value of Il^i 27rn. ^ Details of this renormalization procedure will 
be given in p3[ . 
The result is 



|Det 1/) I 



Here we used the formula 



n 



2 sin 



I det(M 



2 



(31) 
(32) 



n 



x~ \ smTra; 



(33) 



The second step is the diagonalization of ^ . If we represent p using the 
basis (|2^) , it is easy to see that p is the block diagonal matrix. ^ We can obtain 
diagonal representation of ]p by diagonalizing this block diagonal matrix again. 

First we assume that M has no degenerate eigenvalue. If C is a stability 
angle of the orbit, — C and C* are also stability angles of the orbit because M is 
the symplectic matrix . Therefore eigenvalues of M can be classified into four 
types, and M can be transformed into a block diagonal matrix by a symplectic 
transformation. 

/mi \ 

m2 



M : 



ruk I 



1. elliptic: ^ = zba 



2. hyperbolic: C = ±i/3 



cos a 
sin a 







— sma 
cos a 



.-0 



(34) 



(35) 



(36) 



^This rule can also be derived from analytic continuation of a generalized zeta function. 
Let C,a{s) = ^a^", then exp(— ^^(0)) is regarded as a normalized value of ]^^in- In this 
case, a„ = 27rn and Ca{s) = {2n)~''(^{s). is Riemann zeta function.) Therefore Cai^) — 
{2-17)-= {(' (s) - C(s)ln(27r)) and cl(0) = 0.{Here we use ({0) = -1/2 and c'(0) = -Inv^-) 
After all wo obtain exp(— ^^,(0)) = 1. 

■^In general, a real symplectic vector space with a given quadratic form can be decomposed 
into a direct sum of skew orthogonal real symplectic subspaces, and the quadratic form is 
represented as a sum of normal forms on these subspaces. This is known as Williamson's 
theorem ^, ^. 
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3. inverse hyperbolic: C, = ±i(/3 + tt) 



m 



4. loxodromic: ( = ±a ± iP 



m 





cos a 



( cos a 


sin a 
y e~''sina 



-e^ sin a 



(37) 



cos a 



— e ^sina 


e~''cosa y 



(38) 



The matrix elements between a;^,„ and x^/ are non-zero only if ( and ( belong 
to the same block. Therefore it is enough to diagonalize each block separately 
to diagonalize ]p . 

For example, eigenvalues and eigenvectors which belong to a hyperbolic block 

are 

Mx± = e±^a;±. (39) 
Here, the normalization condition is taken as 

[x+,x-] = l, (40) 

and eigenvalues and eigenvectors of D are 

Dx±^n{t) = e±^nX±,n{t), (41) 

(t) = exp{e±^n)V{t)x±{t), (42) 
e±,n = T/3 + 2n7ri. (43) 
We define real vectors a;c,±,„,Xs^±^„ as 

Xc,±,n = -^{x±,n+X±-n), (44) 

v2 

Xs,±,n = —^{x±,n-X±-n) (n = l,2, ...). (45) 

Since x±fi are real vectors, we don't have to re-define them. Matrix elements 
which have non-zero values are 



Jo 

/ dtxl_f_^„l/)Xc,-,n 

Jo 

/ dtxl^j^^^lpXs,-,n 

Jo 



—2n-K, 
2n7r, 



(46) 
(47) 
(48) 
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dtxi^+^^ipxc-,n = -/?, (49) 

I dtx^^Ipx.fi = -(3. (50) 
Jo 

Therefore Ip is represented by the matrix (i„ in the space spanned by a;c.±,n, Xs,±^n'- 



dn = 



In the case of n = 0, 



/ 


\ 2mT 






'2m: 



do 



-fi 2n7r \ 


y 



-p 

-13 



(n > 1). 



The solutions of the eigenvalue equation det((i„ — A/) = are 



A = ±V/32 + (2n7r)2 (n > 1), 
A = ±/3 (n = 0). 



(51) 



(52) 



(53) 
(54) 



The solutions for n > 1 are doubly degenarate. Thus we obtain diagonal repre- 
sentation of Ip for hyperbolic blocks. 

Other blocks can be treated in the same way. First, we define real basis 
from the real parts and the imaginary parts of the basis (|2^) , and represent If) 
using this basis. Then Ip become a block diagonal matrix and each block is 
4 X 4 at most. Second, we diagonalize these blocks, then we obtain diagonal 
representation of Ip . 

Diagonal elements of an elliptic block are 



a + 2n7r (n = 0, ±1, ±2, ...). 
For a inverse hyperbolic block, 

±^JWVi2^^+W (n = 0,l,2,...). 
For a loxodromic block. 



± V^2q^(c^+2^i^ (n = 0,±1,±2...). 



(55) 
(56) 
(57) 



All diagonal elements are doubly degenerate. These results are summarized in 
Table |l|. Q 

^Note that the absolute values of the diagonal elements has no generic meaning because 
we use the non-orthogonal (symplectic) basis. On the other hand, the signs of the diagonal 
elements have generic meaning. 
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Till now, we have assumed that the monodromy matrix M has no degenerate 
eigenvalue. If M has degenerate eigenvalues, we have to deal with these cases 
separately. We will not go into the details of these exceptional cases, but the 
case in which M has eigenvalue 1 is very important. ^ In this case, the stability 
angle of the eigenvalue is 0, which means that If) USA a zero-mode and the 
orbit is not isolated. Therefore we have to calculate the functional determinant 
of Ip in the space where the zero- mode is removed (we denote this as Det Ip ) , 
and integrate with respect to the zero-mode separately. 

We deal with the simplest case in which the eigenvalue 1 is doubly degen- 
erate. (This type is called parabolic.) We can choose the basis Xa {a = 1,2) 
as 

Mxi ^ xi, (58) 
Mx2 = X2+^fXi, (59) 

[xi,X2]^l. (60) 

In this case, the basis ( p5| ) are not complete set. Therefore we define new basis 
as 

Xa,n{t) V{t)e^p[t{2mril - B)]xa (a = 1,2). (61) 

Here, / is the unit matrix, and B is defined as 

Bxi = 0, (62) 
Bx2 = 72^1- (63) 

Note that B is the generator of M{e^Xa — Mxa). D acts on these basis as 

Dxx^n — 2mTixi^n, (64) 
Dx2,n = 2m:ix2M - T^^i^n- (65) 

We define real basis as in (Eq), then Ip is represented in this space as 



/ 2mT \ 

-2mr 

-2m: 7 

\ 2n7r 7 / 



(n > 1), (66) 



*The eigenvalue 1 is always degenerate; The number of other eigenvalues {j^ 1) is even, 
therefore the number of the eigenvalues equal to 1 is even. 
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The functional determinant of this part (except zero-mode) is 



Det'^ = 7]^(27^n)^ (68) 



n=l 



^ 7- (69) 

The integration of the zero-mode must be done separately. 

If ]p is derived from the stationary phase approximation of the time-independent 
Hamiltonian system, it has at least one zero-mode corresponding to time-translation 
symmetry. In this case 7 = —dT/dE and the result of the integration of 
the zero-mode is proportional to the period T: 

K = ^ ■ (70) 

^27rh\§^^\det{M' -I)\ 

Here M denote the monodromy matrix whose parabolic part is removed. The 



factor y 27r?ii^ is canceled by the Fourier transformation (||), then we obtain 
the well-known amplitude factor of the Gutzwiller trace formula: 

T 

(71) 



i?iv^|det(M' -/)|' 

Details of the integration of the zero-mode will be discussed in ||l3|| . 

The third step is the calculation of the Maslov index /Lt. If A = 0, the Maslov 
index is considered to be 0, and we can see relative changes of the Maslov 
index by continuously deforming the operator p . Thus we can calculate the 
Maslov index. The point is that the deformation of Ip corresponds to that 
of V{t). V{t) is a curve on the symplectic group, therefore deformation of p 
means deformation of the curve on the symplectic group. If we continuously 
deform the curve without changing the initial point (origin) and the end point 
(monodromy matrix), the Maslov index of the curve doesn't change because 
the diagonal elements are determined only by the monodromy matrix and don't 
change during the deformation. Therefore we regard such curves are equivalent 
(Fig.^. (We refer this equivalence relation as '-~^.) This means that we can 
regard two operator ]/) (A) and Jp (A ) equivalent if A and A are connected by 
a continuous gauge transformation. 

However, if two curves which have the same terminal points can't be shrunk 
to a point continuously, two curves doesn't necessarily have the same Maslov 
index. Such cases actually happen, because symplectic groups have non-trivial 
topology and the fundamental group of them are ||ll| 

7ri(Sp(2n,i?)) Z. (72) 



In other words, the Maslov index (and the path integral (19)) is the multi- valued 



function of the monodromy matrix. The Maslov index is completely determined 
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by the curve on the symplectic group starting from the origin (we refer this set 
of curves as F), then the Maslov index is uniquely determined on the quotient 
space F/ ^, and this is exactly the definition of the universal covering space; 

F/- =Sp(2n,i?). (73) 

Therefore the Maslov index is regarded as the function on Sp(2n, R), instead of 
Sp(2n, R). This is similar to the construction of Ricmann surface of the complex 
multi- valued function, for example, \ogz. 

A point on Sp(2n, R) is specified by two quantities: a symplectic matrix and 
a winding number. First we investigate the change of the Maslov index which 
comes from the winding number. For example, if we continuously change a to 
a + 27r in elliptic case, the Maslov index changes by 2. Therefore in this case 
the change of the winding number by 1 corresponds to the change of the Maslov 
index by 2 (Fig.|). This rule can be applied to all cases, which can be proven 
as follows. 

Let Al and A2 be curves on the symplectic group. We assume that both 
have the same terminal points and difference of the winding number between 
A\ and A2 is 1. Then deform A\ continuously to A2 through the other curves 
Bl,B2 (Fig.^. The difference of the Maslov index between Al and A2 is the 
same as that of _B1 and B2, because the change of the Maslov index during 
the deformation Al Bl is canceled by that of B2 A2 if the end-point go 
through the same path during the two processes. Therefore the rule that holds 
for some cases holds for all cases, and the change of the winding number by 1 
corresponds to the change of the Maslov index by 2, regardless of the stability. 
This means that (Det iJpY^'^ changes sign by the gauge transformation U that 
changes the winding number by odd number: 

{Det iJp (A^)}i/2 = _{Det iJp {A)}^!"^ (74) 

where AP is defined as 

AV ^-U-^'^^U-^AU. (75) 
at 

A similar result is known in anomaly of SU(2) gauge theory coupled to Weyl 
fermions [p^ . 

Now we know the contribution from the winding number to the Maslov index, 
the next question is the contribution from the monodromy matrix. Since the 
winding number changes the Maslov index by a multiple of 2, the monodromy 
matrix determines whether the index is even or odd. The Maslov index changes 
at the point where ip has zero eigenvalues, and this means the monodromy 
matrix has the eigenvalue 1. (This is the condition for bifurcation.) Therefore 
the Maslov index changes when end-point of V crosses the region defined by 
det(M-/) = (Fig.|). 



11 



In Fig.^, we show eigenvalues of the monodromy matrix in the complex plane. 
The stability changes when eigenvalues collide, and collisions are classified as 
follows: 

1. elliptic hyperbolic 

two eigenvalues collide at 1 

2. elliptic <-^- inverse hyperbolic 
two eigenvalues collide at -1 

3. elliptic(x2) ^ loxodromic 

eigenvalues collide simultaneously at conjugate points of the unite circle 

The Maslov index changes when eigenvalues collide at the point 1, and we can 
see from table^ that the parity of the Maslov index changes when and only when 
the stability changes elliptic ^ hyperbolic. |^ 

A given curve C on the symplectic group can be deformed into a normal form 
UkC using continuous deformation (Fig.^). ^ Uk is the closed curve which have 
the winding number fc, and C is defined not to cross bifurcation regions. Since 
the monodromy matrix can be decomposed into some blocks, C can also be 
decomposed into these normal forms; 



elliptic 



• hyperbolic 



cos at — sin at 
sin at cos at 



e/3t 



(77) 
(78) 



inverse hyperbolic 



2/3(t-l/2) \ 1 

) ^2-'-^^ 

^Strictly speaking, the basis ([25j) are not linearly independent when the monodromy matirix 
has eigenvalue 1. Therefore we must use other basis in this region. 
^Here we define the product of two curves A and B as 

AB(t)-{ (0<t<l/2) 
^^y^l - \ A(l)B{2t -1) (l/2<i<l). 

In this definition, A{l)B(2t — 1) means the product as matrices. 



12 



• loxodromic 

/ e^* cos at 



-e^* sin at \ 

e~^*cosai -e~^*smat 



a0t 



sin at 






e~^* sin at 



s''* cos at 








e ^* cos at y 



(81) 



• parabolic 



1 7t 
1 



(82) 



The elliptic and the inverse hyperbolic type have the Maslov index 1, hyperbolic 
and loxodromic type have 0. Therefore if the monodromy matrix have p elliptic 
blocks and q inverse hyperbolic blocks, the Maslov index of the system is 



H = p + q + 2k 



(83) 



If the monodromy matrix have parabolic blocks, we must add to the 

Maslov index. But the factor corresponding to time translation symmetry is 
canceled by the Fourier transformation. 

We are now able to calculate n-repetition of this or bit. Since C" = UkC UkC . 
can be deformed into UknC we only have to consider n-th power of the nor- 
mal forms. The elliptic type increases the winding number when na goes over 
27r, therefore the Maslov index of this part is 1 + 2[na/2'K\. The square of the 
inverse hyperbolic type becomes hyperbolic type, and the Maslov index increase 
by 1. The cube of it is inverse hyperbolic, the Maslov index increases by 1, and 
so on. Therefore the Maslov index of this part is n. The hyperbolic type, the 
loxodromic type, and parabolic type don't change the Maslov index when we 
raise them to n-th power. Thus the Maslov index of the n-repetition of the orbit 
is 

p 

'na.; 



P 



1 



27r 



qn + 2nk, 



(84) 



where Ui denotes the stability angle of i-th elliptic block. Note that /Lt„ ^ n/xi 
if the orbit has elliptic blocks. 

In conclusion, we have derived the semi-classical trace formula using phase 
space path integral. The point is that quadratic path integral around the peri- 
odic orbit is determined by the symplectic flow around the orbit, and the flow is 
regarded as a curve on symplectic group. We directly diagonalize the quadratic 
path integral, and Maslov indices are defined by the signs of these diagonal el- 
ements. Since they are determined only by the monodromy matrix, two curves 
on symplectic group which have the same terminal points and can be shrunk 
to a point is considered to be equivalent. The set of curves divided by this 
equivalence relation become universal covering space of symplectic group, and 
Maslov indices are uniquely determined in this space. A point in the space is 
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specified by a symplectic matrix and a winding number. The winding number 
changes the Maslov index by 2, and the stabihty of the orbit determines the 
parity of the Maslov index. We also derived the formula for the Maslov index 
of n-repetition of a orbit. If the monodromy matrix of the orbit have elliptic 
eigenvalues, /i„ is not generally equal to n^i, as pointed out by Brack et al.||^ 
The author would like to thank Professor T.Hatsuda for helpful discus- 
sions and encouragements. He would also like to thank Professor M.Sano, Dr. 
S.Sugimoto, Professor T.Kunihiro, Professor H.Kuratsuji, Professor A. G.Magner 
and the members of the Nuclear Theory Group at Kyoto University for valuable 
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stability 


A 


diagonal elements 


pet ^ 1 


index 


elliptic 




a + 2n7r 


4sin''a/2 


odd 


hyperbolic 


e±/3 


±/3,±v//32 + (2n^)2 


4sinh^^/2 


even 


inverse hyperbolic 


_g±/3 


±^/32 + (2n+l)27r2 


4cosh^ /3/2 


odd 


loxodromic 




±^/32 + (a + 2n7r)2 


4(cosh/? — cosa)2 


even 



Table 1: Stability, eigenvalues of the monodromy matirx, diagonal elements of 
Tp , functional determinant ip of the part, and the parity of the Maslov index. 
All diagonal elements are doubly degenerate except ±/3 in hyperbolic case, n 
runs over all integers in elliptic and loxodromic cases. In hyperbolic case, n > 1, 
and in inverse hyperbolic case, n > 0. 




periodic orbit 



Figure 1: The space of the displacement vectors is regarded as fibre bundle, and 
Hamiltonian symplectic flow defines the connection. The structure group of this 
space is Sp(2n, i?). 




Figure 2: Curves on the symplectic group. A and B is equivalent and others are 
not. 
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Figure 3: Flow of diagonal elements of elliptic type when a changes into a-\-2-n. 
The vahics of diagonal elements in left side and the right side are the same, but 
the Maslov indices differ by 2. 





det(M-I)=0 



Figure 5: The Maslov index changes at the point where the end-point of the 
curve crosses the region defined by det(M — /) = 0. 




X 



Figure 6: eigenvalues of a sympectic matrix in the complex plane. 
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Figure 7: The change of the sign of the diagonal element of If) when the stabilty 
changes elliptic to hyperbolic. 




Figure 8: Deformation of a given curve C into a normal form UkC . Uk starts 
from the origin, goes around the hole k times, and returns to the same point. 
C doesn't cross the region defined by det(M — 7) = 0. 
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